Berry phase in Magnetic Superconductors 
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In magnetic systems, electronic bands often acquire nontrivial topological structure characterized 
by gauge flux distribution in momentum(k)-space. It sometimes follows that the phase of the 
wavefunctions cannot be denned uniquely over the whole Brillouin zone. In this Letter we develop 
a theory of superconductivity in the presence of this gauge flux both in two- and three-dimensional 
systems. It is found that the superconducting gap has "nodes" as a function of k where the Fermi 
surface is penetrated by a gauge string. 
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Geometric phases and topology of wavefunctions have 
been the subject of intensive studies jl], |). The quan- 
tized Hall effect (QHE) in a 2D electron system under a 
strong magnetic field can be also interpreted in terms of 
the topological integer called Chern number ||. If the 
system has a gap, the Hall conductivity is proportional 
to the Chern number and is quantized. This nontriv- 
ial topology of Bloch waves is not specific to the QHE, 
but is common and universal in systems with broken 
time-reversal symmetry. Namely, in magnetic materi- 
als with the spin-orbit interaction and/or tilting of spin 
configuration, there occur many band-crossing points 
which are locally described as Weyl fermions and act as 
(anti)monopoles in the momentum-space [0, ||, |[ fjj. The 
existence of (anti)monopoles means that the Bloch wave- 
function cannot be defined in a single gauge choice ||. 
The anomalous Hall effect in ferromagnets represents this 
topological property of the Bloch wavefunctions j|, f?j . 

On the other hand, coexistence of magnetic ordering 
and superconductivity (SC) is found in many materials 
and is a recent significant issue ||. Although detailed 
analysis of each material is still missing in most of the 
cases, it is highly desirable to establish the general feature 
of the SC in systems with broken time- reversal symmetry. 

In the present Letter we shall study the SC made 
out of the Bloch states with nontrivial topology. The 
quantity of central importance is the gauge flux de- 
fined in the momentum space, which is generated by 
the vector potential A(„)(k) = — z(kn|Vk|krt), where 
|kn) is the Bloch wavefunction of the n-th band. Band 
indices are represented by subcripts with parenthesis. 
This A(„)(k) represents an overlap of the two wavefunc- 
tions separated infinitestimaiTy in the k-space, i.e. the 
Berry phase connection. The gauge flux is defined as 
B(„)(k) = Vk x A(„)(k). It is worth noting that the 
monopole density p(„)(k) = ^-Vk ■ B( n )(k) is nonzero, 
and is given by /J( n )(k) = YliQlnSQ* ~ k ;„), where qi n 
is an integer called the strength of the monopole p | . 
The monopoles are located at "diabolical points" JT|, fl0|, 
where bands touch each other. The monopole and an- 
timonopole act as the source and the sink of the gauge 



flux B(k), respectively, as in Fig. |l|(a). In general, the 
up- and down-spins have different band structures; there- 
fore, A, B, and p are spin-dependent, and we shall write 
as A Q , B Q and p a , where a is a spin index. 
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FIG. 1: Distribution of (a) B(k) and (b) B A (k). The open 
and the solid circles represent the monopoles and the anti- 
monopoles, respectively. 

The question we address below is an effect of this gauge 
flux B(k) on the SC properties. In order to character- 
ize topological properties of the gap function A(k), we 
define a vector potential A A (k) = — VklmlnA(k) = 
-V k argA(k), a flux density B A = V k x A A (k), and 
a monopole density p A (k) = ^Vk • B A (k). For the mo- 
ment, we assume that A(k) is nonvanishing in almost 
all over the BZ |l^]; we shall discuss later what happens 
without this assumption. As explained later, distribu- 
tion of B A (k) is confined into strings as in Fig. |l|(b), 
due to bosonic nature of A(k). The main result of the 
present Letter is that in 3D topological structures of the 
gap function A(k) is solely determined from that of the 
wavefunction: p^pfe) = p Q (k) — p@(— k). It reflects the 
fact that A Q/ 3(k) represents a pairing between (k, a) and 
(— k, f3) electrons. Thus if the band structure in the nor- 
mal state is known, we can immediately calculate the 
monopole density p A ^(k), irrespective of details of the 
attractive potential. The monopole density /9 A «(k) tells 
us about zeros and phases of the gap function. 

Henceforth we focus on the SC where only one com- 
ponent of the gap function is nonzero: (i) singlet SC, 

(ii) triplet SC with A u ^ 0, A TT = A u = 0, and 

(iii) triplet SC with A TT + 0, A u = A u = §. 
The case (iii) is appropriate for half-metallic magnetic 
SC. In the singlet (triplet) SC, we have A(k) = A(-k) 
(A(k) = -A(-k)), implying B A (k) = B A (— k), and 
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p A (k) = — p A (— k) [ji3|. We note that the time-reversal 
symmetry breaking is a necessary condition for nontriv- 
ial topology of the gap function. In (iii), broken time- 
reversal symmetry is assumed from the outset. We have 
p A (k) = p|(k)— pi(— k). Thus, inversion symmetry must 
be broken to have nonzero p A . In (i) and in (ii) we have 
p A (k) = p|(k) — pi(— k). Hence in nonmagnetic SC, 
(k) vanishes. In magnetic SC, including both ferro- 
magnets and antiferromagnets, /? A (k) is nonzero. 
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FIG. 2: Monopoles (open circles), antimonopoles (solid cir- 
cles) for A(k), and the Fermi surface. The open squares rep- 
resent point nodes on the Fermi surface. 

Let us consider the distribution of B A (k). B A (k) = 
— Vk x VkargA(k) is nonzero only when A = 0, i.e. 
Re A = ImA = 0, which determines a curve in the 3D 
BZ in general (Fig. (b)). Distribution of B A (k) is like 
a Dirac string, starting from a monopole and terminating 
at an antimonopole. In contrast, B(k) does not share this 
feature; distribution of B(k) spreads over the BZ (Fig. |l| 
(a)). Now we discuss the relation among the Fermi sur- 
face (FS), the (anti)monopoles and the Dirac string. For 
any closed surface S in the BZ, the total flux penetrating 
it is quantized. We take the FS as S for gapless sys- 
tems. When the total flux penetrating one of the FS's 
is nonzero, the gap function cannot be expressed as con- 
tinuous. It is quite contrary to our conventional view 
fl4j| . Because p A (k) = — p A (— k), a FS symmetric with 
respect to k = encloses zero monopoles in total, and 
A(k) can be continuous on this FS (Fig. ||(a)). There are, 
however, two cases where the Dirac strings penetrate the 
FS's. One is the case with broken inversion as well as 
time-reversal symmetry. Then, the FS is not symmetric 
with respect to k = 0, and can enclose monopoles with 
nonvanishing strength (Fig. || (b)). The other is a pair 
of FS's symmetric with respect to k = 0. It is possi- 
ble that one FS encompasses a monopole and the other 
encompasses an antimonopole (Fig. ^(c)). For example, 
if one of the FS encloses a monopole of unit strength, a 
Dirac string starting from this monopole necessarily pen- 
etrates the FS. An intersection of this curve with the FS 
is nothing but a point node. This string can intersect 
the FS more than once, and the number of point nodes 
on this FS should be odd. The nonvanishing total flux 
implies that A(k) cannot be a single continuous function 
on each FS. The FS should be divided into regions, in 
each of which A is continuous. 

Here we remark on "conventional" nodes in anisotropic 



SC. Line nodes appear when ReA and ImA have a com- 
mon real factor /(k) as A(k) = /(k)[Reg(k) + ilmg(k)]. 
Zeros of f(k) determine a surface, whose intersection 
with the FS is a line node. It does not affect B A (k), 
because Vim In A (k) = Vlmlng(k). Hence, only point 
nodes are concluded from the gauge flux argument. Dif- 
ference of nonmagnetic SC from magnetic SC lies only 
in absence of magnetic monopoles; in nonmagnetic SC, 
Dirac strings necessarily form loops, which may cross the 
boundary of the BZ. There are many experimental meth- 
ods to get information on point nodes. In particular, heat 
conductivity measurements tell us about the direction 
of the nodes Jll| . When the positions of the nodes are 
known, care must be taken to assign the phase of A(k) for 
magnetic SC, because our theory asserts that A(k) need 
not be continuous on the FS. Group-theoretical classifi- 
cations of gap functions in Jl(], [ijj do not apply to the 
SC with nontrivial topology, since they assume that gap 
functions are continuous. It would be interesting to clas- 
sify all possible multivalued gap functions. 

To see why topological structure of the wavefunction 
is inherited by the gap function, we note the following. 
If the wavefunction is a continuous function of k, there 
is no monopole. Thus, if there is a monopole, we should 
divide a surface surrounding it into regions, in each of 
which the wavefunction is continuous p8| . This resem- 
bles Dirac monopoles where the vector potential 
cannot be a single continuous function in the whole space. 
This "patch" structure of the wavefunctions remains in 
the gap function via the BCS term in the Hamiltonian. 

To understand this mechanism in detail, we start with 
a 2D model and generalize it to 3D. Similar topological 
arguments as we have developed in 3D are also possible 
in 2D as well. In 2D, we define a Chern number and 
a total vorticity of A(k) as Ch = ^ J BZ e? 2 k£>(k) 2 and 
v = j- J BZ d 2 k£? A (k) z , respectively. We can show that 
the total vorticity v a p of the gap function A a j3 is the 
sum of the Chern numbers for spin a and for spin (3: 
v a p = Ch Q + Chp pfj| . We shall explain how this result 
comes out in the singlet SC on the 2D model proposed 
by Haldane @ as an illustrative example. We take this 
model because it is the minimal model which exhibits 
nontrivial topological structure of wavefunctions with- 
out a uniform magnetic field. It is a tight-binding model 
on the honeycomb lattice. Because there are two sublat- 
tices, the Hamiltonian is written as a 2 x 2 matrix, which 
is conveniently expressed with the Pauli matrices er as 
H(k) = vl + w • er, where v = 2i 2 cos<^> VJ • cos(k ■ bj), 
wi = ti J2j cos(k • u> 2 = tiX^sin(k • a,), and 
W3 = M — 2<2 sin <f> . sin(k • bj), where ai = a(l,0), 

a 2 = a(-i,^), a 3 = -ai - a 2 , b x = a 2 - a 3 , 
b 2 = a 3 ai, b 3 = ai - a 2 . t\, t 2 , M, <p are con- 
stants. The Pauli matrices are not associated with spins, 
but with the sublattice structure. We set M = and 
< 4> < tt. Time-reversal symmetry is broken by a stag- 
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gered flux, while the uniform field is absent. This stag- 
gered flux is not necessarily external; staggered flux can 
be realized by a spontaneous magnetic moment in the 
material itself, through spin-chirality and/or spin-orbit 
coupling [7]]. Let us calculate the Chern number for 
the lower band. By dividing the BZ into two regions V± 
as in Fig. ||, where k± = ±^(i, are two inequiva- 
lent corners of the BZ, we can write down the eigenvector 
of the lower band as v + = 1 / ■ ■ VA " 



A /2|w|(|w|-tu 3 ) \ Wi + iw-2 

for k e V + , and v = i ( ~ m +) W ] 

^/2|w|(|w|+tos) \ W 3 + |w| 

for k £ There remains some freedom in the di- 

vision of the BZ. Because v± is ill-defined only at k!j-, 
respectively, we are free to deform this division as long 
as k° ^ V±. This corresponds to the gauge degree of 
freedom. At k S V + n V-, two choices of the wavefunc- 
tion are different by a phase factor e*^( k 'V+ = v_, i.e., 
A+ - A_ = -V k 0(k), where e 1 ^ = Wl Z iW2 »- Thus, 
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j21| . The nonzero Chern number implies that the wave- 
function cannot be written as a single function for the 
entire BZ. This also affects the definition of field oper- 
ators a(j)k- Second-quantized Hamiltonian is written as 
^ = Ek.tj ■■ c kjHji(k)cki. Field operators a(i) k are writ- 
ten as = J^j v li)jO s -) c j'ki- Let a k ± denote the annihi- 
lation operators for the lower band when k S V± . Then, 
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FIG. 3: Division of the Brillouin zone of the Haldane model 
into two regions V+ and V- . 

We now consider the singlet SC on the lower band. In 
conventional formalism in the BCS theory of SC, it is 
assumed that the field operators are continuous in the 
whole BZ; it is no longer true in the present case. We as- 
sume that the Fermi energy crosses only the lower band. 
The BCS pairing term in the Hamiltonian is 



E 

kev+ 



^ A(k)+a k ^ 



(1) 



Singlet SC implies that A(k) + = A(-k)_ for k € V+, i.e. 
A is an even function of k. The guiding principle to study 
topological structure of A(k) is the gauge invariance, cor- 
responding to the freedom in the division of the BZ. It is 
equivalent to require that the BCS term is continuous at 
the overlap of two regions; we get A(k) + = e 2lc ^ k )A(k)_ 



for k e V+ fl V- . The total vorticity v is 
y-v J dargA(k) 
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Nonzero v implies that A(k) can never be expressed as a 
single continuous function for the entire BZ. To see what 
the gap function looks like, we take an example of an on- 
site attraction: Hjj = —U J2i n iT n il- We assume that 
there are two hole pockets surrounding k^. each. The at- 
tractive potential V (k, k') for the lower band is factorized 
as a product of a function of k and that of k'. We obtain 
thereby A(k)± = C Wl(k g^ 2(k) (k € V±), where C is 
a complex constant. This has two vortices at k = kfj_, 
with vorticity —1 each. The spin Hall conductivity u s xy 
p2| can be calculated as u s xy = 0; it is quantized when 
A(k) ^ on the FS, and a xy can be nonzero in general. 

We can construct simple 3D models starting from the 
Haldane's model. For example, we set wi = t\ J2j cos(k • 
a,-), w 2 = ti J2j sin(k-a J ), w 3 = ~2t 2 cos(fc 2 c) J2j sin ( k ' 
hj) in the Hamiltonian H — w ■ cr, where c is the lattice 
constant in the c-direction. This is a tight-binding model 
on a stacked honeycomb lattice with hopping along ±aj 
(nearest-neighbor) and ±bj ± (0, 0, c); the latter hopping 
acquires ir/2 phase if it is clockwise in the hexagonal pla- 
quette. This model has p = ^ Q (3=± (3S(k — k° + /3^c), 
leading to p A = 2p 7^ 0. In some range of parameters, 
some of the FS's encircle a monopole or antimonopole, 
and the total flux is nonzero for them. 

So far we have dealt with noninteracting systems. One 
may wonder if the topological properties of A(k) dis- 
cussed above are robust against interactions. In general, 
robustness of physical phenomena originated from topol- 
ogy is guaranteed by an existence of integer topological 
numbers, which remain unchanged under an adiabatic 
change of the Hamiltonian. Although this adiabatic prin- 
ciple usually applies to gapful systems, it applies also to 
Fermi liquid with interactions. In Fermi liquid, bare elec- 
trons turn to quasiparticles which are well-defined near 
the FS. (Even with interactions, the FS is topologically 
stable in the momentum space. See Fig. 1 in (23).) The 
BCS theory is founded on this Fermi liquid theory, and 
the gap function A(k) is defined in terms of these quasi- 
particles. The total strength of A(k)-monopoles inside 
the FS is a surface integral of a well-defined function 
-^Vk x V k argA(k) over the FS, and it remains a well- 
defined integer even when we include interactions and 
when A(k) vanishes far from the FS. When we start with 
noninteracting systems, and adiabatically switch on in- 
teractions, this integer remain the same for certain range 
of the strength of the interactions. We note that also the 
total strength of monopoles for the wavefunction inside 
the FS remains an integer in the presence of interactions. 
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Its generalized definition in the presence of interactions Science and Technology of Japan, 
is zfae^Tr fgdS^Gd^G^GdkrG-iGd^G- 1 where 

G(zw,k) is the one-particle Matsubara Green's function, 

and S is the three-dimensional surface surrounding the 
FS in the four-dimensional frequency-momentum space 

It is a well-defined inte- 
pj|. After applying the 



(w,k) (see Eq. (65) in (23)). 
ger, and topologically stable 
Stokes' theorem, the integrand is finite only on the FS, 
and the integral is equal to the number of monopolcs 
inside the FS. In 3D, Fermi liquid picture is valid even 
with interactions, and modification of Fermi velocity vf 
or of quasiparticle residue Z does not change the num- 
ber of monopoles inside the FS. Thus the total monopole 
strength remains unchanged. It changes only when the 
topology of the FS's changes. 

Let us discuss on applicability of our theory to real ma- 
terials. As mentioned above, we should search through 
magnetic SC with (a) equal-spin pairing without inver- 
sion symmetry or (b) opposite-spin pairing. The case (a) 
is likely in ferromagnetic SC, while (b) is likely in anti- 
ferromagnetic SC. Band structure calculations have been 
concentrating on energy eigenvalues on high-symmetry 
directions. Almost no attention has been paid to ex- 
istence of band-touching points, or to phases of wave- 
functions. Thus, we do not have at present sufficient in- 
formation to see which materials are candidates for our 
theory. In order to consider how often nontrivial topolog- 
ical structure appears, we should instead resort to simpli- 
fied models. In @], the authors consider a ferromagnetic 
tight-binding model on the 2D square lattice with ti g or- 
bitals, with physically reasonable values of parameters. 
The gauge flux B(k) z has sharp peaks, where the bands 
approach each other in energy. If we go to the 3D by 
adding an extra dimension along k z , the bands will touch 
each other at some points, because band-touching points 
(diabolical points) has codimension three jioj ^3) . Actu- 
ally in J(|, the diabolical points, i.e. the (anti)monopoles, 
are found in a model for an antiferromagnet. These 
points are topologically stable. Thus nontrivial topology 
of the wavefunctions and the gap functions, as discussed 
in the present Letter, is expected in many materials. 

We can argue exciton condensation in the similar way. 
In 2D, total vorticity of the order parameter is equal to 
the difference of Chern numbers of the bands involved. 
In 3D, the monopole density of the order parameter is the 
difference of those for wavefunctions of the two bands. 

In conclusion, we considered the SC in a band with 
nontrivial topology. It is formulated in terms of a gauge 
flux distribution in the k-space. Although the gap func- 
tion A(k) depends on a detail of attractive potential, 
topological structure of A(k) is determined solely from 
that of the normal-state wavefunction. The effect of dis- 
order is an interesting issue, and is left for future studies. 

The authors thank helpful discussion with R. Shindou 
and M. Sigrist. We acknowledge support by Grant-in- 
Aids from the Ministry of Education, Culture, Sports, 
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